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Loschmidt echo with a non-equilibrium initial state: early time scaling and enhanced 

decoherence 

Victor Mukherjee/'Q Shraddha Sharma, 2 ^ and Amit Dutta 2 <[! 

1 Institut de Physique Theorique, CEA Saclay, F-91191 Gif-sur-Yvette Cedex, France 
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We study the Loschmidt echo (LE) in a central spin model in which a central spin is globally 
coupled to an environment (E) which is subjected to a small and sudden quench at t = so that its 
state at t = + , remains the same as the ground state of the initial environmental Hamiltonian before 
the quench; this leads to a non-equilibrium situation. This state now evolves with two Hamiltonians, 
the final Hamiltonian following the quench and its modified version which incorporates an additional 
term arising due to the coupling of the central spin to the environment. Using a generic short-time 
scaling of the decay rate, we establish that in the early time limit, the rate of decay of the LE (or 
the overlap between two states generated from the initial state evolving through two channels ) 
close to the quantum critical point (QCP) of E is independent of the quenching. We do also study 
the temporal evolution of the LE and establish the presence of a crossover to a situation where the 
quenching becomes irrelevant. In the limit of large quench amplitude the non-equilibrium initial 
condition is found to result in a drastic increase in decoherence at large times, even far away from 
a QCP. These generic results are verified analytically as well as numerically, choosing E to be a 
transverse Ising chain where the transverse field is suddenly quenched. 

PACS numbers: 03.65.Yz, 05.50.+q, 05.70.Jk, 64.70.qj, 64.70.Tg, 75.10.Jm 



The emergence of the classical world from the quan- 
tum world, namely decoherence, or the quantum-classical 
transition through reduction of a pure state to a mixed 
state has been a subject of perpetual interest to the 
physics community^— . The concept of the LE has been 
proposed in connection to this quantum-classical transi- 
tion in quantum chaos to describe the hypersensitivity of 
the time evolution of a system to the perturbation ex- 
perienced by its surrounding^—. The LE is defined as 
follows: if a quantum state \ip) evolves with two Hamil- 
tonians H and H' , respectively, the LE is the measure of 
the overlap given by 

£(£) = |(-0|cxp(i#'f)expH.fft)|V')| 2 . 



In recent years, the temporal evolution of the LE has 
been studied in the vicinity of a QCP. In this context, 
the central spin model (CSM) where a central spin (CS) 
is coupled globally to all the spins of an environment (E) 
which is chosen to be a transverse Ising chain has been 
introduced^. The CS is assumed to be in a pure state 
initially while the spin chain is in the ground state. The 
interaction between the central spin and the environment 
effectively leads to two Hamiltonians which provide two 
channels of time evolution of the environmental ground 
state and lead to a decay in the LE. It has been reported 
that in the limit of weak coupling between the central 
spin and the environment, the LE shows a sharp decay 
close to the QCP of the spin chain and right at the QCP, 
it shows a collapse and revival as a function of time t. 
This collapse and revival of the LE can be taken to be 
an indicator of the proximity to a QCP. It can also be 
shown that the CS makes a transition to a mixed state 
when the LE vanishes. 

The CSM has been generalized to a more general- 



ized environmental Hamiltonianii, to the limit of strong 
coupling^, to the case when the interaction between the 
CS and E is local (i.e., CS coupled to a single spin of the 
E) rather than global^ and also in the realization of a 
Schrodingcr magnet^. In the short-time limit, one finds 
a sharp decay of the LE given by the Gaussian forrn^ 
C(t) ~ exp(— at 2 ), where the decay rate a is expected 
to capture the universality associated with the QCP of 
the E 16 i 18 ~ — . In a recent experimental study with NMR 
quantum simulator, it has been observed that for a fixed 
short time the LE approaches a minima in the vicinity 
of the QCP of an antiferromagnetic Ising chain^. The 
connection between the dynamic LE approach with the 
static fidelity approach has also been established^. 

In this paper, we study the temporal evolution of the 
LE following a sudden quench of the E by a perturbation 
at t = so that the state of E at t = + , is an eigen- 
state of the initial Hamiltonian (Hi) but not of the final 
Hamiltonian (Hp). Our aim here is to study how a sud- 
den quench of the E influence the LE (or the decoherence 
of the CS) both in the early time limit and also in its time 
evolution. Denoting the excited (ground) state of CS by 
|e) (|g)) one can write down generic Hamiltonian of the 
system and environment in the formic 



H(X + 5, g) = H + (X + S\e) (e|) V x + gV g , 



(1) 



Here, H + XV\ is the initial Hamiltonian Hj(X), of the 
E where A is the measure of the deviation from the QCP 
(A = 0) while gV g denote the perturbing term responsi- 
ble for the sudden quench so that the final Hamiltonian 
H_f(A,s) = Hi(X) + gV g ; V\ and V g do not necessarily 
commute with each other. Unless otherwise stated g is 
always taken to be positive. The term <5(|e)(e|)VA is the 
global coupling between the CS and the E where we have 
assumed that the CS couples with the E only when it is 



in the excited state |e) and also the small S limit through- 
out. Initially (t = 0) the CS is assumed to be in a pure 
state, \4>s) = c g \g) + c e \e), with the \c g \ 2 + \c e \ 2 = 1. On 
the other hand, since the E undergoes a sudden quench 
at t = 0, it is in the ground state of the Hamiltonian 
Hi(X) at t = 0+ denoted by \G(X,g = 0)). 

As a consequence of the coupling between the environ- 
ment and the CS, the state \G(X,g = 0)} starts evolv- 
ing following two separate hamiltonians Hp(X + S,g) = 
H + (X + S)V X + gV g and H F (X,g) = H + XV X + gV g 
depending on the state of the CS. At an arbitrary time t, 
we can therefore write the state of the composite system 
as 

|V(t)> = c g \g) ® \4> g (t)) + c e \e) \4> e {t)), (2) 

where \<j> g {t)) = exp(— iHp(X, g)t)\G(X, g = 0)) and 
\<f) e {t)) = exp(-iH F lx + 6,g)t)\G(X,g = 0)). The LE 
at time t is defined as 



c q (x l9l t) = \^ g (t)\Mt))\ 2 

= \{G(X,g = 0)\e lH ^ x ^ t e~ lHF< - x+5 ^ t \G{X,g = 0))K&) 

To contrast with the equilibrium situation (g = 0)^, we 
note that in the present case \G(X,g = 0)} is not an 
cigenstate of H(X,g). Throughout we shall denote the 
LE for the equilibrium case (no quenching) as C. 

Having introduced our model, let us briefly summarize 
our results. In the following, we propose a general scaling 
for decay rate a of the LE in the short time limit close to 
the QCP and show that the scaling is unaffected by per- 
turbation V g . We however show that there is a crossover 
to the equilibrium behavior for a limit of the quenching 
amplitude g when the correction in the LE due to the 
quenching becomes irrelevant. When the quenching is 
relevant, it deeply influences the temporal evolution of 
the LE and shows a faster decay as a function of time 
even when the E is quenched to the QCP and one ob- 
serves a conspicuous absence of the proper revival when 
the environment is chosen to be a transverse Ising chain. 

To arrive at the scaling relation valid in the short time 
limit, we assume the limit S — > and truncate the ex- 
ponentials in Eq. (|3|) up to the order t 2 . One can then 
express the LE in the form C q (X, g, t) w 1 — at 2 ps e~ at , 
where a = S 2 [(V 2 ) — (Va) 2 ] , where (...) implies the ex- 
pectation value in the state \G(X,g = 0))£. This implies 
that the small time limit behavior of the echo is indepen- 
dent of the quench amplitude g. 

Let us note that the operator XV\ is a relevant or 
marginal perturbation that drives the E away from the 
gaplcss QCP thereby generating a gap in the spectrum 2 ^; 
it is then expected to scale in the same way as energy, 
implying XV\ ~ X vz for A ^> L _1 /' y , where v and z are 
the correlation length and dynamical exponents, respec- 
tively, characterizing the QPT in the E driven by A . We 
eventually find the general scaling form 

a ~ ,5 2 A 2 ^- 2 (A > L- 1 /"); a ~ S 2 L 2 ^- 2z (A < 2T 1/!y ),(4) 

where L is the linear dimension of the system. However, 
we note that the above scalings arc valid only as long as 



2/v—2z > d, where d is the spatial dimension. Otherwise 
the contribution from the low energy modes become sub- 
leading and that from the higher energy modes dominate 
resulting in the scaling a ~ L d . In deriving the above 
relation, we have assumed that the CS is coupled to the E 
through the operator Vs which is chosen to be equal to V\ 
that drives the system away from the QCP. Otherwise, 
the critical exponent v§, that determines the scaling of 
correlation length when the system is moved away from 
the QCP (A = 0) due to the perturbation Vs (correlation 
length, £ ~ 5~ l ' 6 ) 1 will replace v in Eq. 0J We note that 
studies for the g = case in the recent past have yielded 
similar results^. 

Let us now reconcile the scaling given in Eq. @ with 
that presented in reference [131 ] when v = z = 1. In 
the short time limit, if one considers only low energy 
modes up to a cut off momentum k (assuming that the 
critical mode k c = 0), we find that the energy gap scales 
k z - L~ z close to the QCP. With A^ A - L~ z , we find 
identical result a ~ A~ 2 L -2z to that given ir>i£ which 
has also been verified for situations with z ^ 1—. Scaling 
relations presented in (Eq. [4j are also in congruence with 
the dimensional analysis noting that 5 ~ A ~ L^ 1 ^ and 
t~X~ vz - L z . 
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FIG. 1: Variation of C q {t) as a function of hi for L = 200, 
S — 0.025, g — 0.2 and t = 1, as obtained numerically from 
Eq. IJTJ. There is a sharp dip around the QCP. 



Let us shift our attention beyond the short-time limit 
to explore the effect of quenching on the LE at an ar- 
bitrary time but for small values of g. Assuming that 
\(dG(X,g))/dg) exists everywhere and using Eq. ((3}, one 
finds 



C q (X,g,t)&£(X,0,t) + g 



d£ q (X,g ,t) l 
dg 



1 9=0 



(5) 



where £(A, 0, t) is the Loschmidt echo in the equilibrium 
situation 5 = when the environment is initially in the 
state |G(A,0)} which is the ground state of Hp(X,g = 
0) = Hi (A). Hence, A c = C q (X,g,t) - C(X,0,t) varies 
linearly in g if dC q {X, g ,t) / dg\ g= o exists and is non-zero. 
Demanding the Loschmidt echo to be dimensionless, we 
find that dC q (X,g, t)/dg\ g= o must have the dimension 
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FIG. 2: Variation of q/(1 — ft) as a function of (1 — h) inside 
the ferromagnetic region for hi — hf = h, g — 0, L = 500, 
t = 0.005 and 8 = 0.005. q/(1 - h) falls as (1 - h)' 1 implies 
a is independent of (1 — h). Inset shows «/|l — h\ varies 
symmetrically on either side of the QCP. 




FIG. 3: The decay rate a as a function of system size L for 
g = 0, hi = h } = 0.999, t = 0.005, 5 = 0.0001 is shown. The 
contribution from the higher energy modes makes a increase 
linearly with L. Inset shows the variation of a/g as a function 
of g for L = 100, t = 0.0001, 8 = 0.001 hi = 0.8, hf = fa + 3. 
a/g falls as <? _1 implying that a is independent of g. 



g^ 1 ~ L 1 /" 9 where f g is the correlation length exponent 
corresponding to the operator V g . We therefore propose 
a very interesting crossover scenario; in case the quench- 
ing parameter g ^$> L~ x ^ Vb , the quenching influences the 
LE as given in Eq. (|3J|. In the other limit g <C L~ 1 / Ug , 
the correction due to the quenching becomes irrelevant 
and one recovers the results for the equilibrium case^ 3 -. 
Eq. (jSJ) also shows that Ac depends on the sign of g, 
i.e., the difference in the LE (compared to the g = 
case) depends non-trivially on the direction of quench 
for small values of \g\. However, in the limit of large \g\ 
(\g\ ^$> L~ l / Vg ), £(A, 0, t) is in general finite far away from 
the QCP; this can be easily seen by expanding C(X,0,t) 
for small <5 and noting that d\G(X, Q))/dX is large only 
near the QCP and small otherwis e) 13 ' 16 . 

On the other hand, in order to investigate the depen- 
dence of C q {\ g, t) on g, we expand G(A, 0)) in Eq. §5§ 
as G(A, 0)) = ^ n a„|S n (A,g)) in terms of the eigen- 



states |S n (A,<7)) of Hp(X,g) (with corresponding prob- 
ability amplitudes a n and eigenenergies e n (X,g)). In- 
crease in g increases the number of oscillatory terms con- 
tributing to the expansion of exp (iHp(X, g)t) \G(X, 0)) = 
X) n exp(ie„(A,0)t)a„|E ra (A,flr)), which finally leads to 
decrease in value of C q (X,g,t) at any finite time. There- 
fore in the limit of large g we can expect A^ to be neg- 
ative at large times irrespective of sign of g, or in other 
words, a sudden quench of an environment coupled to a 
qubit leads to enhancement in decoherence as compared 
to the case with equilibrium initial condition. We note 
that the dependence of C q (X,g,t) and A£ on g as dis- 
cussed above are not necessarily related to a QCP, and 
should be prominent even away from it. This is one of 
the central results of our work and it shows that equi- 
librium initial condition is absolutely necessary for a co- 
herent time evolution in an open quantum system, even 
far away from the QCP. It is worthwhile to note that 
recent studies have shown that non-equilibrium critical 
dynamics (achieved by a slow variation of a parameter) 
of the environment across its QCP also leads to dramatic 
enhancement in decoherence^. 

We elucidate the above generic theories using the ex- 
actly solvable example of a spin 1/2 in an environment of 
transverse Ising spin chain so that the composite Hamil- 
tonian is 



H I (h + 5) = -J2(,°-°. 



j'+i 



he 



6\e)(e 



(G) 



The environmental spin chain undergoes a QPT at h = 
h c = ± l 24 i 25 . At time t = 0, the transverse field is sud- 
denly changed from the initial value hi to the final value 
hf = hi + g so that the initial state is the ground state 
\G(hi,0)) of the Hamiltonian H](hi). We are consider- 
ing the case when the quenching along direction of the 
transverse field that drives the quantum transition (i.e., 
V g = V X ). 

The Loschmidt echo at time t > can be calculated 
exactly in the form 

C q (t) = \(G(hi)\e iir <- h ri t e- iH l h t +s ) t \G(hi))\ 2 = Y[F k 



I \A k e lAkt + B k e~ lSkt - C k e 



iS k t 



k>0 
£> fee -iA fc i|2 



(7) 



fc>0 



where 

A k = cos a k (h i: hf) cos a k (hi, hf + 5) cos a k {h f , h f + S) 
B k = cos otk(hi,hf) sin ak(hi,hf + 5) sin a k (h fl hf + 5), 
C k = sinak(hi,hf) cos a k (hi,hf + 5) sin a k (h fl hf + 5), 
D k = sinakih, hf) sina k (hi, hf + 5) cos a k (hf, hf + S), 
A fc = e k (h f +5)- e k (h f ) and S k = e k (h f + 5) + e k (h f ), 
a k {m,n) = (9 m — 8 n )/2, 9 m = tan^ 1 sin(fc)/(m — cos A;) 

and s k (h) — 2d (h — cosk) + sin 2 k. Wc note that the 
limit g = (i.e., hi = hf), we recover the result of the ref- 
erence [13(, where we have used periodic boundary con- 
ditions such that the decoupled momentum modes k are 
quantized as^ k — 2irp/L, p = 0,1, 2, ...L. 



We shall now analyse Eq. (Jjj) numerically to study the 
proximity to the QCP. Fig. (1) shows that C q {t) plotted 
as a function of hi shows a minimum near hi = h c — g = 
1 — g, thus detecting the presence of a quantum critical 
point even when the environment is suddenly quenched. 

Let us now verify whether the scaling relations of a in 
the early time limit given above holds true in this specific 
case. We present results in Fig. ©; the variation of a 
is plotted as a function of A = 1 — h close to the QCP 
with g — 0; we find the results are in perfect agreement 
with the proposed scaling. In inset of Fig. (2), we show 
that the short-time decay also depends on the symmetry 
of the phase on either side of the QCP—. Moreover, a 
is linear in L which is essentially a contribution from the 
high-energy modes as we show in Fig. ([3]) where we also 
show that the decay rate a does not depend on g for 
.9^0. 

In Fig. (fj|, we plot the behavior of C q (t) as a function 
of time, as obtained numerically using Eq. ([7]). Inter- 
estingly, for nonzero g, L q initially decays for small time 
before rising again and showing small oscillations at large 
times. On the other hand the echo stays almost close to 
unity away from the QCP for g = 0. One therefore finds 
that quenching deeply influences the collapse and revival 
of the LE following a quench. In fact, absence of perfect 
revival (even when the E is quenched to the QCP with 
hf + 5 = 1) and very fast decay of the peak in C estab- 
lishes this conjecture. Moreover, we consider the limit 
5 <C g when sin a k (hf, hf + 5) <C cos a k (hf, hf + 5) and 
one gets 



F, 



\A k e lAkt 



D k e 



-iA fr f |2 



(8) 



Near the critical momentum mode k — 0, A k = £k{hf + 
S) — Sk(hf) ~ 2(hf + 5 — hf) = 25 for g 3> L _1 , eventually 
leading to 

F k « (A k + D k f cos 2 {A k t) + (A k - D k ) 2 sin 2 (A fc t). (9) 

The time at which the maxima occurs can be deter- 
mined by the condition dF k /dt — which yields t n ,k = 
nir/A k ~ <5 _1 (n = 1,2, •••). However, the momentum 
dependent t n ,k results in the irregular oscillatory tempo- 
ral behavior of jC q (t), as shown in Fig. (QJ. In Fig. (|5|), 
we show the time t$ at which the first maxima occurs 
indeed scale with d^ 1 . 

Let us now estimate characteristic value of g = go ~ 
jj-i/vg. f or g > g o ^ we ex p ec t to see the influence of the 
non-equilibrium initial condition. Inspecting, Eq. ([7]). 
we expect a cross-over when the contributions from the 
terms C k and D k become negligible as compared to 
those coming from A k and B k . This is ensured when 
sma k (hi,hf) « cosa k (hi,hf). Using the initial and fi- 
nal parameters to be hi = 1 — g/2, hf = 1 + g/2, we find 
go ~ L , as is expected from our previous analysis be- 
cause v g = 1, in the present case. In the limit of small g 
it is straightforward to verify numerically that the differ- 
ence Ac scales linearly with g as proposed in the generic 
case (see Fig. ©). Interestingly, as expected from our 
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FIG. 4: This figure shows the variation of C q (t) with time for 
L — 100, S — 0.025 and various values of g, hi and hf (as 
marked in the figure). The blue and black curve corresponds 
to g = whereas the red and green curve marks values corre- 
sponding to g — 0.4. The green curve corresponds to crossing 
the QCP Inset: C q {t) with time for L = 100, S = 0.025, 
hi = 0.575 and g — 0.4, so that hf +5 — 1. 
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FIG. 5: Time to plotted as a function of 5 for L = 100, 
hi = 0.95, hf = 1.05. As expected, to falls linearly with 5. 



earlier discussions, even though the sign of dAc/dg is re- 
versed for small \g\ when direction of quench is reversed, 
in contrast in the limit of large \g\ (\g\ ^> go) Ac asymp- 
totically decreases to —1 irrespective of the direction of 
quenching, signifying almost complete loss of coherence 
(see Fig. ©) . 

In conclusion, we have studied the LE both in the early 
time limit and as a function of time following a sudden 
quench of the environment. We have proposed a generic 
early time scaling which is independent of the quench- 
ing and determined entirely by the scaling dimension of 
the operator V$ that couples the CS to the E and the 
dynamical exponent z which is a characteristic of the as- 
sociated QCP of the environment. Moreover, our study 
predicts an interesting crossover behavior. We also show 
that non-equilibrium initial condition can lead to drastic 
increase in decoherence in open quantum systems, even 
away from a QCP. For the specific example in which a 
transverse Ising chain is chosen as the E, we not only 
verify the generic scaling relations but also establish that 
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FIG. 6: Ac decreases (increases) linearly with g for L — 
200, t = 20, 5 = 0.025, ft, = 1 - p/2 and ft/ = 1 + .9/2, 
(hi = 1 + g/2 and ft/ = 1 — g/2) in the limit of small g. 
Inset : Ac asymptotically decreases (to ~ — 1) for large g, 
with L = 200, t = 20, <5 = 0.025, ft* = 1 =F fl/2 and ft/ = 
1 ± <?/2, confirming that for large g, the coherence disappears 
completely irrespective of proximity of a QCP and direction 
of quenching. 



the perfect revival of LE does not occur. However, in 
this example we have chosen V\ = V g = Vs. It would 
be instructive to verify when the CS is coupled to the E 
through a longitudinal field term (<5|e)(e| Yli a f ) or the 
E is quenched by a term g ^ erf . Although the model 
then becomes non-intcgrablc, the exponents v g and us 
will be different from v = 1 (in fact, the critical expo- 
nent v when the QCP is perturbed by a longitudinal field 
is2£ 8/15 ) and are expected to appear in the early time 
scaling and the crossover relation. Further, this work 
can also lead to additional studies on effects of different 
quenching schemes on decoherence and generalization to 
finite temperatures may also yield interesting results. Fi- 
nally, recent experiments on Loschmidt echo using NMR 
quantum simulator— and also other notable experiments 
on non-equilibrium quantum dynamics (for example see 
Ref. [231) have created the possibility of experimental 
verification of our results in near future. 
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